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J/ip production in PHENIX 
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Abstract. Heavy quarkonia production is expected to be sensitive to the formation of a quark gluon 
plasma (QGP). The PHENIX experiment has measured J/ip production at y'sjvjv = 200 GeV in Au+Au 
and Cu+Cu collisions, as well as in reference p+p and d+Au runs. J/tp's were measured both at mid 
(\y\ < 0.35) and forward (1.2 < \y\ < 2.2) rapidity. In this letter, we present the A+A preliminary results 
and compare them to normal cold nuclear matter expectations derived from PHENIX d+Au and p+p 
measurements as well as to theoretical models including various effects (color screening, recombination, 
sequential melting...). 
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PHENIX is one of the experiments located at the Rel- 
ativistic Heavy Ion Collider (RHIC) of Brookhaven Na- 
tional Laboratory. It has the capability of measuring quarko- 
nia through their dilepton decay in four spectrometers: 
two central arms covering the mid-rapidity region of \t]\ < 
0.35 and twice tt/2 in azimuth, and two forward muon 
arms covering the full azimuth and 1.2 < \r]\ < 2.2 in 
pseudorapidity. Electrons are identified in the central arms 
by their Cerenkov rings and by matching the momentum 
of charged particles reconstructed in drift chambers with 
the energy deposited in an electromagnetic calorimeter. 
In the forward arms, muons are selected by an absorber 
and identified by the depth they reach in a succession 
of proportional counters staggered with steel walls. The 
event vertex and centrality are measured by beam-beam 
counters situated at 3 < |ry| < 3.9. For A+A collisions, 
the centrality measurement is further refined by the use 
of two zero degree calorimeters located downstream the 
beams. A detailed description of the PHENIX apparatus 
can be found in yQ. 

This letter is organized as follows: we first set up the 
references by looking at d+Au data (section ^) and com- 
pare the cold nuclear matter effects derived from it to 
A+A data (section^ . In section[31 we compare our results 
to what was observed at the CERN SPS. We then review 
in section 0] the state of the art of the possible explana- 
tions of our observed suppression. Kinematic distributions 
could help distinguishing between models and we present 
our rapidity (sectional and mean squared transverse mo- 
mentum (section |BJ| before to conclude in section [7| 

Note that all A+A J/tp data shown below are PHENIX 
preliminary. 



1 J/ip production in d+Au collisions 

The nuclear modification factor Rab for any A+B col- 
lision type is defined as the J/ip yield observed in these 
collisions, divided by the yield measured in a p+p run and 
scaled by the average number of nucleon-nucleon collisions 
(Ncoii) extracted from a Glauber model: 
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Departure from Rab = 1 implies some nuclear modi- 
fications of the produced J/^'s. The PHENIX collabora- 
tion has measured the J/tjj nuclear modification factor for 
d+Au collisions 0. Its centrality dependence is shown on 
figure n for our three rapidity ranges. 

In d+Au collisions at 200 GeV, J/ip's in our three ra- 
pidity ranges probe the following momentum fractions x of 
gluons in the gold nucleus (neglecting the emitted gluon): 
0.05 to 0.14 (upper panel, negative rapidity, gold-going 
side), 0.011 to 0.022 (midrapidity) and 0.0014 to 0.0047 (lower 
panel, positive rapidity, deuteron-going side). The differ- 
ence observed between forward (lower panel) and back- 
ward (upper panel) yields indicate that shadowing and/or 
anti-shadowing are at play. These effects are modifica- 
tions of the parton distribution functions in nuclei with 
respect to the usual functions in free nucleons. In partic- 
ular, one expect a depletion (shadowing) at low x due to 
overlapping and recombining partons (gluons in the case 
of J WO- The strength of gluon shadowing is not heavily 
constrained by theory and models predictions differ 
by a factor of three. 

Comparison with Vogt's theoretical predictions 6 shows 
that a modest amount of absorption (0 to 3 mb), added 
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Fig. 1. J/ip nuclear modification factor Rja as a function of 
centrality (given here by the number of collisions N co u) for 
backward (top, y = —1.7) mid (middle, y = 0) and forward 
(bottom, y = 1.8) rapidities. The histograms are the distribu- 
tions of N co u for each of our four centrality classes. Theoretical 
curves are from Vogt and take EKS shadowing and 0, 1, 2 
or 3 mb normal absorption cross section (from top solid to 
bottom dot-dashed). 



to a modest amount of shadowing such as the one pro- 
vided by the Eskola-Kolhinen-Salgado (EKS [3]) scheme, 
can describe both the rapidity and centrality dependen- 
cies, as shown on figure ^ 



2 Normal nuclear matter effects in A+A 



To extrapolate the effects seen in d+Au to A+A col- 
lisions, one has to rely on a model. Together with our 
Au+Au (~ 0.2 nb _1 ) and Cu+Cu (~ 3 nb _1 ) measure- 
ments |3, Au+Au predictions from such a cold nuclear 
matter model jH] are presented on figureEl as a function of 
centrality (displayed here as the number of participants). 
They account for shadowing (following the EKS prescrip- 
tion) and nuclear absorption (1 mb for solid lines and 3 mb 
for dashed lines). Predictions are shown for two rapidity 
values (y = and y = 2) corresponding to our two ra- 
pidity regions. The difference between the 1 mb and the 
3 mb curves give an idea of our current uncertainty on the 
cold nuclear matter effects affecting the J/ip production 
in A+A collisions. 

However, in both rapidity ranges, our most central 
Au+Au measurements depart from the predictions made 
with the stronger nuclear effects (3 mb), suggesting that 
other suppression mechanisms are involved, namely, an 
anomalous suppression. 
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Fig. 2. Nuclear modification factor Raa as a function of cen- 
trality (given here by the number of participants N par t) for 
Cu+Cu y = (open magenta circles) Cu+Cu y = 1.7 (full 
red circles) Au+Au y = (open green squares) and Au+Au 
y = 1.7 (full gold squares). Error bars are statistical, brackets 
are systematics, and box are global (common to all the points). 
Theoretical curves are Au+Au cold nuclear effects only predic- 
tions from Vogt 0, solid and dashed lines being for 1 mb and 
3 mb normal nuclear absorption cross sections, respectively. 



3 Lower energy anomalous suppression 

Such an anomalous suppression was early predicted by 
Matsui and Satz as a signature of the QGP |5j and later 
observed by the CERN NA50 experiment ^01 a t lower 
energy (^/snn — 17.3 GeV). Various models explain the 
NA50 data and the extrapolation of three of them to 
RHIC energies are displayed on figure 01 together with our 
Au+Au data. 

In (solid lines), J ftp's are absorbed by comoving 
particles (of undetermined partonic/hadronic nature). In 
|12j , the authors describe the dynamical interplay between 
suppression and regeneration of J ftp's in a QGP. The sup- 
pression mechanism is dominant for NA50 energies and is 
the only one presented here (see figure CD for the full pre- 
diction). In JSj (dot-dashed lines), the authors use a QGP 
statistical charm coalescence model. All three models fail 
to reproduce our data, overestimating the measured sup- 
pression. Other models such as percolation J3| a l so over- 
predicts the suppression, suggesting that new mechanisms 
are at play at RHIC energy. 



4 Alternate explanations for suppression 

For now, at least three classes of models exist that can 
roughly accommodate the amount of anomalous suppres- 
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Capella et al. hep-ph/0505032 
Grandchamp et al. hep-ph/0306077 (direct) 
Kostyuk et al. hep-ph/0305277 (QGP) 
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Fig. 3. Nuclear modification factor Raa as already presented 
on figure|2l Only Au+Au data are shown (y — open symbols, 
y = 1.7 full symbols). Theoretical curves are predictions from 
models that describe the anomalous suppression seen by the 
NA50 experiment at lower energies |11H2I13") . 



sion seen in our most central preliminary results. The ac- 
curacy of the present data, as well as the uncertainty on 
the cold nuclear effects, do not allow to favor one or the 
other. We review them in the three following subsections. 



4.1 Detailed transport 

One paper [T^l, simulating J ftp transport in a hydrody- 
namical model, actually predicts an amount of suppression 
that matches our most central data. It is shown on figure^ 
where the authors have added cold nuclear matter effects 
(nuclear absorption only, 1 or 3 mb) with respect to the 
published paper. The suppression they obtain is not large 
because they authorize J ftp's with sufficient momenta to 
freely stream out of the plasma. 



4.2 Sequential melting 
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Fig. 4. Nuclear modification factor Raa as already presented 
on figure|2l Only Au+Au data are shown (y — open symbols, 
y = 1.7 full symbols). Theoretical curves are from a J ftp trans- 
port model in an hydrodynamical quark gluon plasma |15j ■ The 
dashed (solid) lines are with (without) anomalous suppression. 
The up (down) curves are with 1 mb (3 mb) normal nuclear 
absorption. 



4.3 Recombination 

At RHIC energies, multiple cc pairs are produced (10 
to 20 in central collisions according to 53)- Quark mo- 
bility in a deconfined medium could allow uncorrelated 
charm quarks to recombine when the QGP fireball freezes, 
forming quarkonia, and raising their yield with centrality. 
A balance between suppression and enhancement could 
lead to the intermediate suppression we observe. Figure 
shows a collection of predictions from various recombina- 
tion or coalescence models, from |12U18U19II20| . Unfortu- 
nately, the lack of knowledge of the yield and distribu- 
tions of charm quarks initially produced as well as of the 
recombination mechanism, make these predictions hardly 
predictive. A good way to search for hints of recombi- 
nation is then to look at its impact on the distributions 
of kincmatical quantities, such as rapidity and transverse 
momenta. 



An important fraction (30 to 40 %) of J ftp's comes from 
feed-down decays of charmonia excited states (tp' , Xc)- 
This particular point is taken care of in most of the previ- 
ous approaches. But recent lattice computations indicate 
that J ftp's could melt at a much higher temperature than 
the one that was originally thought. One possible hypoth- 
esis, defended in is that only the excited states melt, 
leaving all the initially produced J ftp's untouched. This 
30 to 40 % suppression could also match our data. 



5 J /if? rapidity profile 

According to [50], the rapidity distribution of the recom- 
bined J ftp's should be modified with respect to the di- 
rectly produced J ftp's. In particular, it should emphasize 
the region of phase space where more c quarks are pro- 
duced, because recombination probability goes quadrati- 
cally with the c quark densities. Recombined J ftp's should 
then have a narrower rapidity distribution. Figure[S]shows 
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Fig. 5. Nuclear modification factor Raa as already presented 
on figure|2l Only Au+Au data are shown (y — open symbols, 
y = 1.7 full symbols). Theoretical curves come from various 
recombination models (solid 12., dotted (18) . crossed (19) and 
dashed P7|). 



our measured rapidity spectra for three Au+Au and four 
Cu+Cu centrality classes, as well as for p+p. 

Within our limited statistical accuracy, we do not see 
any modification of the rapidity profile and cannot draw 
any conclusion towards or against recombination from this. 



6 J/ip mean transverse momentum 

Recombination modifies transverse momentum as well. 
Figure shows our measured mean squared transverse 
momentum (jpj,) for various systems and centrality inter- 
vals, as a function of N co u. The shaded yellow bands are 
predictions from |2()j . corresponding to either J/ip's from 
recombination (lower band) or to directly produced J/i/>'s 
(upper band). To properly predict the modified pt spec- 
tra, one first needs to quantify the normal pt broadening 
coming from cold nuclear effects (Cronin effect). This ef- 
fect was clearly seen by PHENIX at forward rapidity, 
by comparing (p^) in p+p and d+Au collisions. In [23j the 
authors use these measurements to quantify the Cronin 
effect. However, at midrapidity, because of the very poor 
statistics of the p+p sample, it is difficult to quantify it. 
Thus, the theoretical prediction can only be safely com- 
pared to the forward rapidity case. 

The solid red line is another parametrization of the 
Cronin effect from |21j. using this simple parametrization: 



(Pt)aa 



(Pt)pp + P<7$((Pt)) x L = 2-51 + 0.32 x L (2) 



where (jij) is in units of (GeV/c) 2 and L is, in fermis, 
the average thickness of nuclear matter seen by a J/ip 
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Fig. 6. J/ip rapidity spectra for three Au+Au centrality in- 
tervals (gold squares), four Cu+Cu centrality intervals (red 
circles) and p+p collisions (black triangles). 
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Fig. 7. Mean squared transverse momentum versus N co u. 
Black stars at N co u — 1 are for p+p collisions, for mid (up- 
per) and forward (lower) rapidity. Triangles are for d+Au, for 
mid (open) and forward/backward (full) rapidity. Squares are 
for Au+Au for mid (open) and forward (full) rapidity. Yellow 
shaded bands are from 1201 and stands from direct (lower band) 
or fully recombined J ftp's. The red line is a parametrization of 
Cronin effect derived from d+Au data in |21|. the dotted lines 
on each side reflecting the associated errors. 
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in the collision system and the centrality bin considered. 
The factor pa5((p^)) = 0.32 is determined by the d+Au 
forward data and stands for the nuclear density p, times 
the elastic gluon-nucleon scattering cross section er, times 
the average px kick given at each scattering 8{{p^)). The 
dotted lines on each side of the main curve are the er- 
rors derived from the p+p and d+Au uncertainties. This 
parametrization provided a good description of lower en- 
ergy experiments |21| and matches our data when it can 
be safely applied, namely at forward rapidity. 

Thus, the raise of (p^) at forward rapidity in Au+Au 
collisions seems explainable by the Cronin effect only. A 
better p + p baseline is needed to interpret the rather cen- 
trality independent (p^) we observe at midrapidity . 
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7 Conclusions 

PHENIX preliminary results show a substantial amount 
of J/ip suppression, growing with centrality. At the high- 
est energy density probed so far, this suppression is as 
large as a factor of three. However, its interpretation is 
not easy. First, the amount of normal, cold nuclear mat- 
ter suppression is not precisely known and demands more 
d+A data. Nevertheless, the observed suppression in A+A 
seems to exceed the maximum suppression authorized by 
cold nuclear matter, and thus be partly anomalous. 

It also seems weaker than the suppression derived from 
models that fit lower energy data. Three classes of models 
that raise the number of surviving J/t/j's were presented. 
They all suppose the formation of a QGP. To distinguish 
between them, a better precision on data, and in particular 
on the kinematical distributions, is required. 
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